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ABSTRACT 

An investigation of the primary shock front and Mach Y-stem 
systems was conducted utilizing decigram explosive charges. Distance 
and arrival time data of the primaty shock front was correlated with 
that of high yield explosions. <A good correlation would indicate the 
feasability of conducting laboratory scale tests to obtain information 
on high yield explosions without the expenditure of time and money 
involved in large scale field tests. The shock front system was also 
investigated at the time of first formation of the Mach Y-estem. 
The critical angle of incidence of the primary shock front for the 
formation of the Y-stem was compared to the theoretical value. 
Theoretical calculations of the yield of the explosion was also comp- 
ared to the actual yield. 

Basie data were obtained from the explosions by photographing 
the shadow of the shock front system utilizing a Polaroid camera, 
a microflash unit and a time delay generator. 

Correlation of the overpressure was excellent for small to mod- 
erete distances from the point of explosion, as was the yield comper- 
ison. The critical angle of incidence comparison indicated a marked 


difference between experimental end theoretical velues. 
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SYMBOLS AND ABBREVIATIONS 


Speed of sound under preveiling ambient 
conditicns 


Ratio of primary shock front radius to 
arrival time 


Calculeted energy of explosion 

Horizontal distance from ground zero (point 
on blast table directly under charge) to 
triple point 


Horizontal distance from ground zero to 
Mach Y-stem 


Geometric factor to convert measured param- 
eters from plane of grid screen to plane of 
charge 

Transmission factor for the speed of sound 
Transmission factor for density 

Height of charge above blast table 

Height of triple point above blast table 


Radius of primary shock front 


Delay time or arrival time of shock front 
at an observed radius 


Arrivel time of sound at an observed radius 
Difference between t, and t 

Charge weight 

Reference charge weight 

Penteaerythritol tetranitrate 
Trinitrotoluene 

Reduced energy parameter 

Reduced radius parameter 

Reduced time paremeter 


Yield factor 
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1. Introduction 

When decigram charges of chemical explosives are detonated, a 
shock front originates at the point of detonetion and expands spher- 
ically outward. At a critical engle of incidence of the primary shock 
front on a reflecting surface, a Mach Y-stem begins to form. This 
critical angle is dependent upon the Mach number of the incident shock 
front. The Y-stem is formed when the primary shock front combines with 
the reflected shock front. As the radius of the shock front continues 
to increase, the Mach Y-stem grows end travels out from the point of 
detonation approximately perpendiculer to the reflecting surface. This 
phenomenon can be photographed at different times and distances from 
the point of detonation and thus, by compiling the data obtained from 
many successive detonations, a study of the primery shock front send the 
Mach Y-stem can be achievede From this study there should result a 
correlation between small chemical explosions and large point source 
(nuclear) explosionse A photograph of a well formed Y-stem is shown 
in Fige 1. A detailed drawing of the aenniona shock front system is 


shown in Fige Ze 
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Schematic of Reflected Shock Front System 
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Figure 2 





20 Experimental Setup 


Ae 


The experimental setup had a two fold purpose. First, it 


contained the necessary electrical circurity to initiate the explosive 


cha rgee 


Second, it provided a meens to photograph the resulting shock 


front eat a preset, accurately mesured time deley. The following 


components were used: 


CY) 
(2) 


(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 


1391-B Time Dealy Generator - General Radio Company 
§24B Electronic Counter with 526B Time Internal Plug-in 
Unit - Hewlett Packard Compeny 

504 Cathode Ray Oscilloscope - Textronix Incorporated 
1530-A Microflesh Unit - General Radio Company 

110B Poleroid Cemera with Tripod 

Oscilloscope Camere = Dumont 

Cathetone ter 

Two Phototubes - RCA type 868 


120 Volt Variac - General Radio Company 


(10) Blest Table 


(11) Scotch-brite Reflective Grid Screen 


(12) Two Lab-jacks 


(13) Warning Buzzer 


(14) Test Light 


Be 


To accomplish the above, the apparatus was positioned and 


connected electrically as shown schematically in Fige 3-6 


To better understand the necessity of this particuler setup a 


brief description of the experimental procedure is given. 


The basic electronic components were photocell #1, the time delay 


generator, end the microflash unit. Photocell #2, the CRO and the 
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electronic counter were only necessary to determine the true time delay 
between the explosion and the light flash. 

The charge was mounted above the blest table and its alignment 
checked with the cathetometer. The shutter of both the object camera 
and the CRO camera were opened manually in a darkened roome The charge 
was then detonated by closing the switch of the variec which was con- 
nected to an AC wall outlet. The blast of the explosion was sensed by 
photocell #1 which sent the start impulse to the time delay generator, 
to the start circuit of the electronic counter, and to the trigger 
circuit of the CRO. The impulse to the time delay generator was de- 
layed a prescribed number of microseconds and then continued to the 
microflash unit to initiate the light flash. This flash of light cast 
a shadow of the shock front on the reflective grid screen» The shadow 
was recorded by the object camera. The light from the microflash unit 
was sensed by photocell #2. The impulse from this photocell stopped 
the electronic counter and initiated a vertical blip on the CRO. The 
true time interval between the explcesion and the photographing of the 
shock front (when the light flashed) was displayed on the electronic 
counter. The time interval was also recorded on the CRO camera film. 
The camera shutters were closed efter the cherge detonetion and light 
flasho 

With this operation in mind, the assembly of the apparatus can 
now be discussed. 

C. Charge, Reflective Screen am Blest Table: See Fige 4. 

The cherge of PETN was mounted at the desired height above the 
blast table by attaching the wires of the charge to the two leads 


from the firing circuit. These leads were passed through a 1/4 inch 
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2e Camera 
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Figure 4a 
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Figure 4b 





copper tube which was at first suspended vertically from a horizontal 
cross arms the cross arm was well above the ares in which the shock 
fronts would be photographed. This method of suspension, however, was 
discarded in fevor of a diagonal suspension due to the slight distortion 
amd possible attenuation of the shock front as it passed vertically up 
the copper tube. Measurement of the primary shock front redius was 
facilitated by having a pure shock front in the vertical plane directly 
above the point of explosion. With the diegonal suspension and the 
flexibility provided by the leads, the charge position could ensily be 
adjusted for exact alignment . 

A cathetometer was used to ensure exact alignment of the charge 
in two dimensions. The desired position of ths charge was at a given 
height directly above a hole centered in the blast table and in line 
with the microflash unit and the center vertical grid line of the 
reflective screen. This position was accuretely determined before the 
charge was mounted by the use of a metal centimeter rule placed up- 
right with orm edge directly over the center of the blest table hole. 
The cathetometer was positioned and leveled, but offset horizontally 
from the center line. This offset was necessary due to interference 
with the microflash unit and camera which were on or near the center 
line. ‘the height of the cathetometer was adjusted to the exact desired 
height of the charge by aligning the crosshairs of the cathetometer with 
the appropriate centimeter mark on the metal rule. The offset of the 
cathetometer also required that a bench-mark be accurately inscribed on 
the grid screen in line with the vertical crosshair of the cathetometer 
and the vertical edge of the metal rule. Therefore, when mounting the 


charge it merely had to be positioned so as to coincide with the bench- 





mark on the screen and the crosshairs of the cathetometer. The cathe- 
tometer wes also utilized to check the level of the blast table by 
moving the metal centimeter rule to the front, to the back end to both 
sides of the blast table and sighting on the desired centimeter mark 
of the rule with the cathetometer. The table was level when the de- 
sired centimeter mark on the rule coincided with the horizontel 
crosshair in the cathetometer end the horizontel grid on the reflect- 
ive screen. The reflective screen was leveled when initially mounted 
against the wall. The blast table (80 by 85 cm.) was supported on its 
sides by two "angle iron" beams which were in turn attached by movable 
clamps to the blast table freme. This freme was given added rigidity 
by the use of numerous cross and diagonal support rods. Two lab jacks 
were placed under the center portion of the blest table. These jacks 
made fine adjustments in height and leveling of the tabie easier as 
well as to prevent the blast teble from buckling under the force of 
the explosive blast. This entire spparatus was placed on a solid 
leborstory teble. 

To reflect as much light as possible to the object camera, a 
screen of highly reflective Minnesota-Manufacturing Company Scotch- 
lite tape was constructed. The tape, in two foot wide strips, was 
attached to a sheet of 4 x 8 foot plywood. To furnish a grid on the 
screen for distance measurement purposes, black thread was attached 
horizontelly and vertically at exact intervals of ten centimeters. 

The completed reflective screen was then mounted and leveled flush 
against the wall. 

D. Firing Circuit. (See Fig. 5) 

The variac was connected to a standard wall outlet (115 volt, 


AC). The veriac rheostat was set at 60 volts to ensure adequate 
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current through the charge wire. For ssfety purposes an extra switch, 
in addition to that on the variec, was added to the circuit. This 
switch wes closed just prior to opening the camera shutters on each 
firing and opened immediately after the room lights were turned on 
after the explosione The variac switch wes used to finally complete the 
circuit. After the object camera shutter was opened, the veriac switch 
was closed and immediately opened. The resulting instantaneous current 
was sufficient to detonate the cherge. The leads from th safety 
switch were passed up the side of the blest table freme end through 

the copper tube to the charge. 

Eo Basic Electronic Circuit. (See Fig. 6) 

Photocell #1 was supported directly below the hole in the center 
of the blast table. To provide continuity in the blast table, so that 
the hole would not interfere with the reflected shock front, a lucite 
rod was placed up through the hole until it was flush with the top 
surface of the blast table. The rod also prevented unexploded PETN 
particles, resulting from low a@der detonations, from emtering the 
photocell. The rod extended down into the photocell terminating 
slightly above the phototube. With this arrangement the photocell 
could sense the explosion of the charge without interfering with the 
shock front. 

The delay generator was positioned well behind the microflash unit 
and in close proximity to the CRO and electronic counter. The lead 
from photocell #1 was conrected to the "PFR Drive" connection of the 
delay generator. The following positions and settings of all electronic 
components were not in accordance with the preliminary set-up proce- 


dure found in the applicable instruction manual. 


il 





The positions and settings for the delay generator were: 

(1) triggering kbvel knob - orange dot in "1030" position and 

white errow on “positive going - AC". 

(2) sweep trigger switch on "direct". 

The output from the delay generator came from the “deleyed sync. 
out" connection am went to the power unit of ths microflash. This 
unit was best positioned on top of the delsey generator. The power unit 
was removed from its case so that its input from the delay generator 
cmild be connected directly across the microphone receptacle, S02, 
(refer to wiring diagram of microflash manuel). The voltage to fire 
the flash was passed from the power unit to the lamp unit using the 
high voltsge cable provided with the microflash unit. Extreme caution 
was exercised in mking the above connections am in manipulating the 
microflash unit since the power ami lamp units contained voltages up 
to 1500 ami 8500 volts respectively. 

The signal from photocell #1 was also connected to the “start 
input" receptical of the electronic counter. This signal started the 
counter counting in tenths of microsecoms. The important positions 
and settings for the counter were: 

(1) Sep-Com switch to "SEP" 

(2) Start and stop trigger slope to "" 

(3) Stert amd stop trigger level voltage to "8x3" 

(4) Display time knob to the 12 o'clock position 

Reliable operation of the electronic counter could only be ex- 
pected after a werm up tim (master switch to "on") of a minimum of 
two hours. The counter was left conected to the wall outlet with the 
mester switch in the down position at all tims when not in use. This 


provided continuous heating of the counter. 
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Lastly the signal from photocell #1 also went to the “external 
trigger in" connection of the CRO. For operation of the CRO the posi- 
tions and settings were: 

(1) trigger: slope switch to "4", coupling switch to "AC", source 

switch to "ext.". 

(2) sweeptime/om variable switch to cslibrated position (full 

clockwise )o 

(3) vertical sensitivity to 20 v/em. 

(4) AC-DC-Gnd switch to "DC". 

The flash of the microflash was sensed by photocell #2 which was 
mounted on a movable ring-stand adjacent to the lamp unit of the micro-= 
flashe A small hole was cut through the aluminum foil at one edge of 
the lamp face through which the photocell could see the flash. This 
hole was shielded so that the light from it would not directly illumi- 
nets the reflective screen. It was determined that for best results in 
stopping the electronic counter, the end of the photocell should be 
placed one inch from the lamp face. 

The signal from photocell #2 was supplied to the stop input connec- 
tion of the electronic counter to stop the time measurement end leave 
the exact elapsed time between explosion end light flash displayed on the 
counter in tenths of microsecondse The only error that could be present 
in this tim display was the error internal to the electronic counter 
itself and the error introduced in the tims it would take for the out- 
put signal from the photocell to build up to the threshold voltage 
required to actuate the electronic counter. The former was eliminated 
by ensuring through calibration that the error was within the design 


Limits of + one microsscond. The latter error was determined to be one 





to two microseconds. This was accomplished by supplying a signal from the 
photocell simultaneously to the trigger circuit and verticel deflection 
circuit of the CRO. The photographs of the CRO traces showed clearly a 
rise time of the photocell signal of one to two microseconds. 

The true time interval as displayed on the electronic counter was, 
therefore, accurate to #three microseconds. This time, however, was 
consistently 35-40 microseconds greater than the time set on the delay 
generator dial. This error was either in the delay generator or the 
microflash unit. By connecting the CRO across the input and output of the 
delay generator, its error was determined to be within the design 
limits of ors percent of the dial setting. The delay between the input 
signal to the microflash power unit and the flash of the lamp varied 
from 35-40 microseconds. This was determined by supplying the output 
from the delay generator simultaneously to the power unit and to the 
trigger input of the CRO. A photocell was positioned to sense the flash 
and to provide its signal to the vertical sweep of the CRO. The time 
delay was read froma photograph of the resulting CRO trace. 

The signal from photocell #2 was also supplied to the vertical 
sweep connection of the CRO. When this impulse was received by the CRO 
it was registered on the cathode-sray tube. This trace which was initi-g 
ated as a horizontal line at the time of explosion and distorted verti- 
cally at the time of light flesh was recorded on the film of the CRO 
cameraco The time was read from this film to an accuracy of about two 
percent. It was therefore merely used as a cross check for the elece 
tronic counter. 

Fo Cemera and Lamp Unit: (See Figo 4) 

The position of the camera was varied depending on the area of the 


shock front that was of primary interest. For the initial pictures 
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of the primary shock front the camera was placed directly beside or 
above the lemp unit. This was satisfectory for the preliminsry quali- 
tative work since it was desired that the picture include the entire 
shock front for all time delays. However, to mike more accurate meas- 
urements of the shock front radius the camera should be moved as close 
to the screen and charge as possible without interfering with the shock 
front or casting shadows on the screen. To eleminate shadows, the 
camera was moved to one side of the center line. This position proved 
unsatisfactory due to the interference of the bright reflection of the 
explosion. In an attempt to move this reflection where it would not 
snterfere, which was to the point on the screen directly behind the 
charge, the camera wes positioned on the centeriine but above the hori- 
zontal line from the light source to the charge. This was accomplished 
by suspending the camera below « tripod. The legs of th tripod were 
extended so that they did not cast shadows on the area of interest. The 
height of the camera was varied so as to keep it as low as possible which 
would keep the explosion reflection low but yet not cast its own shadow 
onthe area where the shock front was expected to appear. It should be 
remembered that the position of the shock from on the film is completely 
independent of the position of the camera since the camera is photograph- 
ing the shadow of tne shock front on the screen end not the shock front 
itself in the plane of the explosion. It should be pointed out that the 
field of view of the camera need nob’ include the center point of the 
blast on the screen due to the mthod used to measure the shock front 
radius. This will be discussed in a later section. 

The optimum position of the camera then, for the first set of 


photographs in which the primary shock front was of interest, was 
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horizontally on the center line but above the charge and approximately 
three feet from the screen. 

It was found that for the masurement of the Y-stem data the camera 
should be moved to one side slightly to give a better view of the Y-stem 
aréa and yet not move the blast reflection into the area of interest. 
Close-up photography of the Y-stem formation was attempted with the use 
of attachable close-up camera lenses. For these close-up photographs 
the camera was moved into approximately 24 inches from the screen. 
sons usable pictures were obtaired, but the most well defined shock 
from images were obtained at the cameras minimum distance of three feet 
without the use of close-up lenses. 

In order to meke the microflash essentially a point source of light, 
the face of the lamp unit wes covered by a square piscs of aluminum 
foil with a one centimeter hole cut in its center. <A point source light 
ensured that the shadow of the shock front was narrow and well defined. 

The lemp unit was placed horizontelly ina direct line with the 
charge and the center grid of the reflective screen and vertically so 
that the light aperture was the exact height of the charge. This posi- 
tion facilitated the measuremsnt of the primary shock front radiuso 
The lamp unit was lowered so that the aperture was level with the blast 
table when the formation of the Mach Y-stem was investigated. This was 
necessary so thet the shadow of the Y-stem would be projected onto the 
reflective screen rather than onto the blast table. 

The distance of the lamp unit from the screen was varied to 
determine the optimum distance for the sharpest shock front. It was 
First thought that this would be accomplished if the light beams were 
parallel, thet is, the lamp placed at an infinite distance. This 


distance was simulated by placing the lamp unit the maximum distance 
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away possible, 18 feet, with negative results. The next position tried 
wes 300 cme This position wes still not satisfactory end the 200 cr 
position was tried. Some good shock waves were photographed at this 
distance but the pictures were too small for accurate messurernt se 

It was decided to move the entire set-up closer to the screen. The 
blast table was moved in flush with the screen which made the charge 
41.5 cm from the screen. Tle light was placed eat 109 cm and the camera 
was placed at three cm. At this light distance the shock wave appeared 
fuzzy even though the grid screen was in sharp focus. Therefore, the 
light was moved back to 165 cm. Good definition of the shock wave was 


obtained at this distance. 
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30 Preparation of charges. 

In preparing the charges for the experiment, a steel die was 
obtained from the Explosive: Laboratory which would compress about 
250 mg of pentaerythritol tetranitrate, PETN, around an approximately 
95 mg detonator. The diameter of the bowl of the die, see Fig. 7, 
was 675 cm which resulted ina spherical charge of .225 cc. The load- 
ing density was 1.5 grams /cce 

A charge was prepared by cutting e 15 cm length of nichrome fuze 
wire in which an overhand knot was tied and adjusted so as to be in the 
center of the wire. This length of wire weighed 25 mg. The knot was 
pulled tight enough so the loop was one mm in diameter. This knot was 
necessary to give a base for the application of the rest of the deto- 
Nator o 

Lead styphnate was used as the first layer of the detonator be- 
cause of its high sensitivity to heat. The lead styphnate was mixed 
with a small amount of Duco cement and acctone. After thorough mixing 
and when sane of the acetone had dried to obtain the proper consistency, 
a small drop of the mixture was applied to the knotted fuze wire. When 
the acetone completely dried, the lead styphnate bead, with the glue 
binder, adhered to the fuze wire and was ready for application of lead 
azide. However, if insufficient glue was used, the resulting styphnate 
bead would flake away or breek apart in subsequent handling. Too much 
glue led to low order explosions, therefore, the minimum amount of 
glue was used to adequately bimi the styphnate. The total amount of 
glue and styphnate used avereged five milligrams. 

A layer of shock sensitive lead azide was applied around the 


styphnate bead using the same mixing procedure as before, only, due to 
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the cohesive nature of the azide itself, very little glue was required. 
Approximately 25 mg of lead azide was used which made the total deto- 
nator weight approximately 55 mg. Since the fuze wire weighed 25 mg, 
the explosive portion of the detonator weighed 30 mg, only 8% of the 
entire explosive cherge. 

It should be noted thet in the preparetion of the detonators, 
actual weighing of each component used was not necessary. Of the 
completed detonator, the weight of wire (25 mg) was fixed. The re- 
maining portion of the detonator was mostly lead azide. So little 
syt phnate was used that its weight varied only within one milligramo 
With the weight of these two components constant, the azide was applied 
until the total detonator weight exceeded 55 mg. This final desired 
weight was used due to the high incidence of low order explosions from 
previous charges when detonators weighed from 35 to 45 mgo With 
practice this final weight of 55 mg was obtained with very small errar. 

The PETN charge itself was made by capsulating the detonator ine 
side the 250 mg of PETN and compressing it in the die. The PETN was 
donated by the Trojan Powder Company of Allentown, Pennsylvania. 

Since it was stored under water, it required being dried for about 

three hours with occasional stirring. If heated excessively it was 
noted that the PETN turned light yellow. The use of this PETN led 

again to low order explosions. 

Once dried, the powder was stored in a dessicator and only left 
exposed to the atmosphere long enough to prepare a series of ten to 
twenty charges at a time. Occasionly, low order explosions occurred 


which were due, in part, to the hygroscopic property of PETN. 





The method of capsulating the detonator in the PETN was as follows; 

ac 285 mg of PETN was accurately weighed out. The use of the 
Mettler type balence greatly decreased the time required to make these 
weighings. 

be The upper and lower sections of the die were placed together 
and about one-third of the PETN was poured into the bowl of the die. 
A small funnel was used to keep spillage at a minimume The PETN was 
tamped down lightly in the bottom of the bowl. (A SAFETY FACE MASK 
WAS ALWAYS WORN DURING THIS OPERATION) 

co The upper section of the die was removed. No PETN was spilled 
if the preceding operation was done carefully. The fuze wire was 
crimped on both sides of the detonator so that a right angle was formed. 
The wire ends were inserted through the small holes located on each side 
of the bowl pedestal and pulled, from the bottom, until the detonator 
rested on the partially filled bowl. It was mcessary to ensure thet 
the wires laid in the grooves provided so the upper section of the die 
was able to slide down snugly over the pedestal. 

ad. When the upper and lower half were again placed together, the 
remaining PETW was poured into the bowl and the plunger inserted. The 
die and plunger were now placed in the press and the press screwed down 
until hand tight. If too much pressure was applied the plunger struck 
the metal bowl and became belled-in. This required either remachining 
or reforming, using a small ball bearing. Also, excessive pressure 
caused the plunger to sever the fuze wire. The die was left in the 
press approximately one minute. 

e. Upon releasing the screw pressure and removing the plunger and 


upper section of ths die, the fuze wires were pulled out of the pedestal 
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base holes. Due to the adhesion between charge and bowl, the charge was 
then released by raising slightly the small brass plunger which protruded 
through the bowl pedestal. This was done by turning a small hex-screw 

on the bottom of the die. Then the charge came out freely. 

fo. The completed charge was eccurately weighed to ensure that the 
total weight was 250 mg plus the weight of the detonator. The weight 
was usually accurate within 5 mg. If the charge weight exceeded the 
desired weight a spatula was used to scrape the excess PETN from the 
charge. If the weight was only slightly less than that desired, the 
charge was used but its actual weight was recorded. 

This procedure wrked quite satisfactorily but was perfected only 
after considerable effort was spent on (1) trying to break the adhesion 
between charge and bowl without losing portions of the charge and there- 
by ruining its spherical geometry and (2) trying to eliminate the spillage 
which resulted when pressure was applied to the plunger. 

In an attempt to eliminate the first problem, liquid graphite and 
liquid mold release was sprayed on the bowl and allowed to dry. The 
applicetion of both failed to release the charge satisfactorily. 
Finally, by enlarging the diameter of the brass plunger seated in the 
bowl, the complete charge came out as described in peragraph (e) above. 

The problem of spillage was never eliminated, although the volums 
of the cavity in the plunger was enlarged am the knife edges of the 
plunger were made sharper, in hopes that the entire quantity of PETN 
would be compressed around ths detonator with the minimum spillage. 
Spillege, however, could not be eliminated with the existing die and, 
therefore, to obtain consistently 250 mg of PETN the problem was circum- 


vented by adding excess PETN. By trial end error it was determined 





thet 35 mg excess was necessary to consistently obtain 250 mg of PETN 
in a completed cherge. The spillage hed to be cleaned from both 


sections of the die prior to preparing snother charge. 
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4, Experimentel Procedure 
Prior to firing a series of charges, e1l electronic equipment was 
energized for a minimum of two hours. After the prescribed warm-up 
period, the electronic counter was checked for timing accuracy by 
switching the Function Selector to "100 KC CHECK" end then to "10 MC 
CHECK". With the Function Selector in each position, the Frequency 
Unit Switch wes cycled through the five positions of frequencye At 
each frequency the counter column corresponding to that frequency 
should register 1.000. If not, the counter required further werm-up. 
The oscilloscope was checked to ensure thet the "Sweep Time /Cm" 
wes at the proper position to ensure meximum use of the entire scope 
face for the prescribed delsy time. This was desirable to give the 
greetest accuracy in reading the delsy time from the CRO picture. 
The "DC BAL" was calibreted as per CRO 504 instruction manual. In 
addition the horizontel position was checked to ensure that the trace 
started exactly at the left edge of the grid. This was accomplished 
by viewing the CRO as a test light wes manually passed across photocell 
#1, thereby producing a test trece. The test light procedure also 
provided a check thet there were no open circuits since it also initiated 
the microflash. It mist be recognized, however, that the trace time did 
not correspond to the delay set into the deley generator. This error 
was caused by the long duration of exposure of the test light to photo- 
cell #1. At times, the menuel sweep test did not initiete the microflesh 
evon though the microflash circuit was completo. In these cases, the 
sensitivity to the microflash power source required adjusting or the 
internal circuitry of the delay generator required cycling. This was 


done by turning the RESET switch to "SWEEP" end then back to "DELAY" 
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severel times. This problem, however, was not encountered when photo- 
cell #1 was initieted with the explosive charge. 

Since the blast table, microflesh, and cethetometer could have 
been moved inadvertently since the last series of firings, the pre- 
scribed distances of each unit from the reflective grid screen were 
checked as well as the level of the blest table. 

The oscilloscope camera and the object camera were loaded with 
new film, if needed, and the focus and lens aperture of the object 
camera were adjusted. 

After the above procedure was completed, the charge wes attached 
to the firing circuit leads, after first ensuring that: 

(1) the variec switch was off, 

(2) the safety switch wes open, 

(3) A FACE MASK WAS WORN. 

The charge was then connected to the leads by twisting the cherge wire 
ends around the leads sufficiently to provide electrical continuity 
end so the charge would not fell off during positioning. 

To accurately position the cherge, it must be moved about in the 
plene of the charge holder until the crossheirs of the cathetometer, 
the charge itself, and the bench merk on the reflective grid screen 
were all in line. 

Once the charge was in position, the object camera was cocked, 
ensuring the exposure time setting was on "B", ami the safety switch 
closed. The variac was set at 60 volts. The room was darkened and 
the apertures of both cameras were opened. A warning buzzer was actuated 
for five seconds and then the variac switch was’cycled to "on" and 


immediately back to the "off" position. This action detonated the charge. 
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Immediately efter the explosion, the apertures >2f both cameras 
were closed, the safety switch opened and the room lights turned on. 
An estimetion was mde et this tim as to the order of the explosion. 
If the exploding charge seemed to spark, as in Fig. 8 or if the sound 
was not of high intensity, or if smell particles of PETN were found 
on the blest table, these symptoms indiceted e low order, incomplete 
explosione However, if the sound of explosion was loud end sharp, 
and a sphere of quickly disappearing exploding ges was observed, and 
there was no evidence of un>xploded PETN, this was a good indication 
of a high order explosion. 

The electronic counter resding and eny of the above symptoms, if 
pertinent, were recorded. The time interval as shown on the CRO film 
was also recorded. In this manner a set of three tines were obteined: 
the deley time set, the counter time, and the CRO time. 

Usually, the series of charges were detonated in rapid order. 
For subsequent firings, it was only necessary to reset the counter, 
check the cameras, dial a new deley generator setting, and attach 
another charge. 

At the completion of the series, all equipment was de-energized, 
except the electronic counter which was always left plugged in but 
with power off, thereby allowing the counter heeter to meintain its 


internal circuitry at a constant temperature. 
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5. Data 

The following procedure wes used to obtain usable data from the 
object camera photographs. 

After the polaroid film wes treated with preservative and allowed 
to dry it wes mounted on a traveling microscope obtained from the USNPGS 
Physics Department. However, magnification wes too great to enable the 
shock front to be distinguished from the background. Since these were 
the only microscopes available and there were no interchangable eye 
pieces or objective lenses available with a lower magnification, a 
negative lense was attached by small pieces of wax directly below the 
objective lens, thereby decreasing the magnification of the film end 
rendering the shock front visible. 

The following data was measured from the photographs and re- 


corded; the primary shock front radius, R,; the height of the triple 


p3 
point, Hips the horizontal distance of the triple point from ground 
zero, Dips and the horizontal distance of the Y-stem from ground zero, 
Die (Table 1) 

The best accuracy in the determination of this data was obtained 
when messurements were made to give the proportion of the distance 
between the two adjacent grid lines between which the shock front was 
observed. The number of whole grid squares from the blast wave to the 
projected point of explosion were then added to this distance obtained 
from the measured ratio. As an exemple, assume the shock front radius 
was approximately 48 cm. The only masurements taken from the photo- 
graphs were the distance from the fourth grid line to the shock front 
and the distance from the fifth grid line to the shock front. This 


retio of distances was multiplied by 10 om and then added to 40 cm to 


oi 





give the actual radius of the shadow of the shock front. This method 

of measurement was the most accurate since the camera could be positioned 
close to the grid screen so thet the distance to be measured included 

the maximum aree of the film. In fact, this method was mandatory for 
accuracy in the lons delay time, large radii explosions due to the 
foreshortening of the grid squares on the film. This error was con- 
siderable for radii larger than 30 cm. The vertical foreshortening of 
the grid square being measured was virtually elirinated by adjusting 

the height of the camera to approximately that height et which the 

blast wave was expected to appear. 

The distances obtained thusly for Ro» Hep» Dé ps Dy, were the 
distances for the shadow of the shock front on the grid screene These 
distances were divided by the appropriate geometric factor to obtain the 
magnitude of the parameters in the actual shock front. tlhe geometric 
factor is dependent only on the distances between source light, charge 
and grid screen. The calculetion of this factor is shown for the last 


physical set-up used. A side-on view of the physical set-up is shown; 


Source Light 


Grid Screen 


~~ 165.8 om Blast Table ‘* 
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From this diagram: 





1ES.8f 124.2 


therefore, the geometric factor, f= ioe = {2S 





This factor applies to Rp as well as to Dy,, and Dy when the 
light was positioned level with the chargee A slightly different 
form of the ratio was used for Htp, since it was measured from the 
blest table to the Mach Y-stem. 

When the source light was placed level with the blast table the 
following diagram and ratio were applicable for Rye A simple ratio 


similiar to that shown above was used for Hyp, Dy, and Dy. 
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Figure 10 
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TABLE I 
Omissions in the table ere due to those parameters not being of 
main interest and, therefore, not in the field of view of the object 
camera. Shots prior to #'42 were for h=5 cm am those following 


#142 were for h=8 cm. 


30 


Shot Noo = _Ry t a? Dep De 
1 19.26 245 04 7032 20017 21.38 
3 19.84 242 08 5044 19.96 20671 
4 19.62 244 63 6056 19637 20017 
6 19.55 242 02 = 18054 20023 
7 19071 24108 6.89 19637 20044 
10 17085 1942 6029 16.81 17.80 
12 22.680 29169 -- = = 
13 22075 29509 8677 21.52 22.670 
14 11.47 84.9 -< a 1054 
15 16.40 16507 -- on 14.52 
ay 24 260 340 00 10522 23016 25028 
18 29070 448 02 15023 26 063 32046 
19 33 060 54504 16076 30055 34 096 
20 39020 646 08 23.00 32.063 39.092 
21 43.00 744 6 22 020 37088 41.96 
23 14.70 13305 4.06 13045 13099 
25 14.70 133 02 == oo 13045 
27 14.00 128 60 -~ c= -- 
39 21.3 26500 on 19.90 20098 
40 19025 185.0 os = a 
105 oe 102.0 2280 11.58 12038 
113 = 88.0 2018 10.29 10.64 
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Shot Noo 


120 
123 
124 
125 
127 
133 
135 
137 
138 
142 
144 
145 
146 
151 
152 
154 
155 
156 
L157 
158 
159 
162 
163 
Wests 
168 
L169 


170 


10576 
10.82 
11.69 
12.26 
12.40 
14.78 
15.45 
16022 


16072 


cow 


18.94 
19eli 
16.56 
14.87 
15.73 
17043 
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oe 


19631 


15.90 


8009 

8304 

88.0 

89.0 

95.20 
14421 
151.5 
164.0 
17303 
279695 
3351.0 
581.0 
25160 
211.0 
171.0 
147.1 
156 20 
190.0 
252 00 
285.0 
52500 
22000 


162.0 
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1.14 
1023 


0.59 


1.86 
BeaG 
3.09 
467 
2026 
2015 
1.86 
0.65 
Iai 
2904 
2671 


4024 


9214 

981 
10.29 
10.08 
10.68 
10.29 
10.24 
10.22 
11.22 
12.91 
13.04 
13.298 
14.69 
21008 
25007 
26038 
17044 
17032 
14 64 
13039 
13.79 
15085 
19.08 
20.04 
22002 
17.81 
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9095 
10.29 
10.29 
10.66 
10.71 
10.41 
1049 
11.52 
1272 
13004 
14.58 
15.03 
21.09 
25016 
26 673 
17.51 
17.52 
15.03 
1363 
13.79 
15085 
19.08 
20.84 
22002 
17.81 
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6.0 Calculations 

Several correlations and calculations were determined from the 
data obtained from the photographs. The most important of these were; 

(1) the Mach number, M, from which the overpressure was calculated, 

(2) the angle of incidence of the shock front when the Mach Y-stem 

initially formed, 

(3) the scaled distances and scaled time, 

(4) and the determination of the average yield of the PETN charge. 

To eliminate the possibility of humsn error in the msny repetitive 
caleulstions described below, a Computer Deta Corporation Model 1604 
computer was utilized. The progrem used is described in Appendix I. 

Ae Mach Number 

The Mach number is the retio of the shock wave velocity to the 


velocity of sound. Thus, a reletionship between distance, time end 


ANG 
M = “Kh = r/dd 
d kp /d £, 


or M dd = teal. 
but f by. = J Ref, 


therefore d é = of Ke Ao. 


by substracting, we have 


di-Jé= te (/- x) 
or M = 


M was obteined, 





/ 
om CA of (at) 
ba 


where pee = ae —dd 
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The quantity lah die, is simply the slope of the curve 
resulting from a plot of at versus Rpe To obtain this slope, the equation 
of the best curve through theat and Rp data first had to be determined. 
(See Fig. 11) The computer gave a set of polynomials from first order 
to sixth order which gave the best fit with experimental deta. By 
comparing the curves of these polynomials with the experimental curve 
it was determined that the fourth order curve most closely agreed with 
the experimental dete. The equation of the fourth order curve used is: 

at = -1.488 + 30.204 R, ~ 09657 R2 + 01889 RB — 00015 RS (2) 

By differentiation of equation (2), the equation of the slope 
of this curve was determined and utilized by the computer to give the 
slope at any point on the curve. This equation is: 

slope = 30.204 — 2(.9657)Ry + 3( .01889 )Re a 4( 00015 )R- (3) 

By substituting the celculsted slope for various values of Rp 
in equation (1), the Mach number wes calculsted. The overpressure was 
then calculated using equation (4). 

p= 2(m*-1) 2 (4) 
where Fo is the atmospheric pressure. 

B. Angle of incidence 

To calculate the angle of incidence, @, it was necessary to 
determine the time at which the Y-stem began to form. This was done 
by using the date from a series of close up photographs with small 
delay time settings. When the data from these shots was incorporated 
with previous data, a smooth curve of height of Y-stem, Hips versus 
time resulted. (See Fig. 12) By extrapolating this curve to zero 
height, the time of Mach Y-stem formation was determined. For a height 
of charge, h, equal to five cm the time was 60 microseconds. For a 


height of charge equal to eight cm the time was 135 microseconds. 
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Using these times, the angle of incidence wes calculated by two 
methods. The first method incorporated the interpolation between 
experimental points of Rp and t to determine Rp for 60 microseconds 
and 135 microseconds. This interpolation gave values of 8.87 cm and 
14.3 cm, respectively, for the radius of the primary shock frmt at the 
time of Y-stem formation. Then, 

= aos lh, 
For the second method, the horizontal distance from ground zero to the 
Mach Y-sten, Dy, wes obtained from the Dy versus t plot at the appro- 
priate times. (See Fig. 13) For t#60 microseconds, Dab 8 cm and for 


t©135 microseconds, D,©12.5 cm. Then, 


© B= tat & 


Shock 


Ab 
Front 3 = 5 or Bem 


By mthod (1), B is 5%.9° end 56° for 5 cm end 8 om, respectively. 
By method (2), B is 60.29 and £7.49 for 5 cm and 8 cm, respectively. 

C. Sceled Distance end Sceled Time 

Scaled distance and scaled time were required for correlation of 
the experimental date with data for any other reference explosion. The 
effects of en explosion vary with the intensity of the blast and with 
the density of the transmitting medium. Since the effect is spherical 
or volumetric a cube root sceling law was used. 

The sceling law for distance is: 


sealed distance f (actual distance) (5) 
cA» 
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a iat {4 mtranusmission factor for densitys (/p Yond we yield factor = 
4) 

Since the experimentel and reference data was obtained at sealevel, 
P=, and fq is equal to unity. Wo is the yield of the reference 
explosion and is equal to one ton of TNT at 70° F. The average experi- 
mental temperature was 71°F. This difference was considered negligible 
am 70°F is used throughout the calculations. To determine tho average 


yield per charge, , in TNT equivalents, the following basic data wags 


g 


compiled: 

Average weight of PETN per charge=249 mg 

Average weight of lead azide and lead styphnate per charges =30 mg 
The equivalents of the above explosives to TNT are 173%, 40%, and 39%, 
respectively. (Reference 1) Using a weizhted average for the TNT 
equivaleats for lead azide and lead styphnate of 39.2%, the total 


equivalents of PETN/charge was calculated. 


Y= 249 + so (44) =25 5.8 ws PETN/chovge 


or in units of tons of TNT, 


W= 225 53 x 193 wy 99 210” eine THT Jehavge 
4AS3.6* 2000 


7 
and, ww = (4.99% /07?)F = ,00722 


The time scaling law is; 


fs x f,x (actual time) (6 ) 
td 


scaled time — 


whers fq and W are defined above. f, is the transmission factor for 
the speed of sound and is equal o(%,). This ratio is proportional to 
the square root of the temperature retio, therefore, fg also is equal 


to oneo 





D. Average Yield 

The scaling law for distance shown above (Eq. 5) can bs used in 
an inverse way to determine the average yield of explosive charge. 
The average yield was calculated to ascertain if the experimental data 
did, in fact, correspond to the known yield of 255.8 mg PETN equivalent. 
The equation in the proper form is; 

WW = Ww 4 (2 senna F 
= Wo “dl \Scaled distance 

where all symbols are defines as aboveoe 

The distance ratio used in Eq. 7 was determined by two methods: 

(1) Visual weighing of the experimental data using plots of b 
ratio versus the log of distance for the experimental and reference 
explosions © 

(2) Analytic determination using the distance - arrival time 
ratio, b, to select a scaled distance. 

Bach method was calculated assuming a point source explosione 
This assumption is justified since the shock front radius data exceeds 
ten times the charge diameter after which the effect of the displacement 
of the atmosphere by the explosive charge mass is considered negligible. 

Method (1), as described in reference 2, is the correlation of 
b ratio with actual distance and scaled distence. In this method 
calculated b ratios were plotted on the ordinsete and the log of distance 
in centimeters was plotted on the abscissa. A similiar plot, using 
the b ratio data from a reference explosion of one ton of TNT (Table 
1l of Reference 2) and the corresponding scaled distance in feet was 
mede. The two plots were then placed on top of each other, so that 
lines of equal b ratio remained superimposed. The two graphs were 


moved horizontally until the curves agreed as closely as possible. 
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The actual distance in centimeters and the scaled distance in feet 
corresponding to this point of intersection were then read from the 
respective abscissae The two distances are corresponding distances 
representing equal b ratio for the two explosions. Incorporating 
these two distances in equation (7), resulted in a average yield of 
275 mg of PETN. Figure 14 shows these two plots on one graph super- 
imposed to give the best fit. 

Method (2), also outlined in reference 1, used the raw data 
directly to obtain the b ratio. the b ratio is tebulated in reference 
2, in feet per millisecond and thus required conversion of Rp from 
centimeters to feet and t from microseconds to milliseconds. The 
simple but tedious calculations were also programmed end printed by the 
computer. By entering Table 11 with the b ratio, a scaled distance in 
feet was obtaimd. This distance was the distance corresponding to 
R end resulted in a surprisingly increasing trend in yield when 
calculated by equation (7). (See Fig. 15) However, when yield is 
extropolated to zero distance, a value of 212 milligrams is obtained. 
This trend is discussed in the next section. 

Ee Data Correlation 

To correlate the data obtained with that from larger scale 
explosions, certain calculations were necessarye The basis for these 
calculations is the energy of explosion, aA, of PETN. BA has been 
calculated for PETN to be 7269 x 10° ergs per graeme For the size 
charges used, 255.8 mg, a@ is 1.91 x 10/” ergs. Using the 1/3 power 
scaling law a reduced energy paremeter, & , in centimeters, was cal- 


a 
oY, = (4/,,) — 26.09 cm 


culated. 
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A reduced radius perameter, )\, and reduced time parameter, 7, were 


also calculated, 


ee Ae and T = at 


OX as 
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7e Conclusions 

The shock front photographs were excellent for the pur pose of 
recording the primary shock front end the Mach Y-stemdatae Jt wes 
observed that at low time delay settings, the reflected front could 
not be seen due to the expaniing geses of explosion. At larger time 
delays, the reflected front was well beyond these gases and was readily 
observed. Occasionally s portion of the secondsry shock front was 
observed. This shock front follows a rarefaction wave into the PETN 
gases. It starts with zero strength and grows as it moves invard 
through these gases. It is swept outward in space until the expansion 
of the high explosive products is nearly exhausted, and then it implodes 
on the origin am is reflected outwerd to the comtact surface. At 
the time it strikes the contact surface between the PEIN gases and the 
air, the PETN gases are still more dense and much cooler then the air 
immediately outside, (this is generally true st eny time). As a 
consequence, the shock in passing through the surfnce sets up an 
inward rerefaction wavee This rarefaction, like the initial rarefaction, 
is followed by a third shock front moving inward. This shock, like 
the second implodes on the origin, reflects snd moves out in the wake 
of the previous shocks. The succession of shocks continues in this 
manner until the energy in the explosion product gases is dissipated. 
(see Ref. 3) The third and succeeding shock fronts were not observed. 
Figure 16 is included from Reference 3 to show the position of the 
shock fronts described above as a function of time. 

The delay tims recorded from the electronic counter were, as 
stated in Section 2, accurate to fthree microseconds. Therefore, 


the basic data R, and t, used in arriving at the conclusions 
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described below is considered a true representation of the actual 
experimental phenomena observed. 

Ae Correlation with high yield explosions 

One avenue of imberest pursued was to determine how close the 
experimental data from the decigram charges of PETN correlated with 
the date from ore ton of TNT. If the correlation was close, the 
feasibility of using laboratory scale explosions instead of tons of 
high energy explosives at a remote test site to simulate the latter 
would be proved. 

Of primary concern in judging the effectiveness of an explosion 
is the peak overpressure, AP, at a given distance from the explosion. 
Figure 17 is a plot of the log of the overpressure versus the log 
of the reduced redius,N. ‘The experimental data obteined for PEIN 
is compared to the theoretically calculated data for TNI obtained 
from Brode (Reference 3). As can be seen, the correlation is quite 
good for the low values of A. At large reduced radii, the overpressure 
is greater than that indicated by the point source, ideal gas curve. 
This ideal gas curve defines the mximum overpressure for a given radius, 
since for an explosion in real air the molecular dissociation and 
ionization of the explosive gases at high temperatures reduces the 
blast efficiency. The discrepency is due to experimental error at 
large radii and can be explained thusly: The light fran the micro- 
flash cast a shadow of the point on the surface of the spherical 
shock front that is tangent to its pathe At low radii this point is 
in the vertical plane through the point of explosion. But as the 
shock front radius increases this tangent point moves out of the plane 


of the charge and toward the light. By simple geometric construction 
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it can be seen that this will give an apparently larger radius. 

Figures 18 and 19 of reduced and scaled paramters, respectively, 
also show good correlation with TNT data from Reference‘’s 2 and Se 

Be. Comparison of known yield with calculated yield 

The result of the yield calculations, although not consistent 
for each method used, gave a calculated yield of the same order of 
magnitude, 275 mg end 212 mg. The per cent error from the known value 
of 255.8 mg of the two methods was 7.5% and 13.2% respectively. 

The first method of calculation as discussed in Section 6 is 
considered more direct and by superpositioning of the two curves 
(Figo 14) many experimental uncertainties are minimized. For this 
reason, more credence is given to the value of 275 mgo 

The linear increase of yield obtained in the second method 
resulted from experimental error which is a function of increasing 
Rp» This method of yield determination should give a horizontal 
straight line at a constant yield of 255.8 mg. The positive slope 
of this line is cmtrary to physical fact. The random error due to 
small variations in charge size and positioning, etc. will flucuate 
som what, but should give an average around a constant yield. By 
extrapolating the curve obtained to zero radius it can be assumed 
that the experimental error, which is a function of Rp , has been 
minimized and therefore the value of 212 mg of PETN is obtained. 

The sources of the experimental error were considered to be a 
combination of several factors, but after considering each factor, 
all except one were discarded. 

The most probable error considered was tim error. However, 


after a careful investigation of the error across all electronic 


41 








components with the CRO, as described in Section 2, the unaccountable 
error was determined to be three microseconds. This source of error 
was therefore eliminated. 

The buckling of the blast table was considered. Although some 
buckling was possible, it was considered extremely small due to the 
leb jack supports located directly under the area of the blast table 
from which the front was reflected. However, any buckling would be 
due to the absorption of energy from the shock front system and hence 
tend to decrease the calewlated yield. This factor, therefore, has 
the reverse effect on yield as thet observed and was eliminated. 

Another source of error considered was the variation in atmos- 
pheric conditions and charge weight. But, the atmospheric conditions 
were reletively constant. The variation of charge weight, although 
small, could possibly have caused an increase in yield of the weights 
increased with an increase in delay times. However, the selection of 
charges was completely random and any error due to this variation 
would have cancelled out. 

The one possible error which did effect yield with an increase 
in Rp was the geomtric factor error as discussed above. This error 
resulted in Ry's greater than actual. When inserted into Equation 7, 
it has the effect of increasing the numerator and decreasing the 
denominator sime the scaled distance varies inversely with the b ratio. 
This term is then cubed which greatly magnifies this error. Therefore 
this error is a function of Rp. By extrapolation to zero Rp, at oe 
assummed that the effect of this error is minimized. 

The values of 275 and 212 mg of PETN yield are considered good 


due to the inherent random nature of experimental explosive researcho 
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C. Angle of incidence for Mach Y-stem formation 

Another avenue of interest was verification of the angle of 
incidence for the formation of Mach Y-stem. From reference 2, the 
angle of incidence is 39° for a renge of Mach numbers from 1.6 to co. 
The calculations from the experimental data (h, Rp, t) indicated a 
rather large discrepency (average angle of 58.1°). However, this angle 
was calculated by two different methods at two different heights of 
charge, h, amd the results were fairly consistent. The geometric 
factor error does not enter into these determinsetions due to the small 
times and radii involved. Further investigation of this incidence 
angle was attempted oat a charge height of three centimeters. 

Clear photographs of the shock front system could not be obtained 
at this height due tothe obliteration of the shock front by the 
explosive product gases at the radius of interest. 

The magnitude of the geometric factor error described above was 
checked by calculeting the corrected Ry using the formla; 

Rp(corrected )=X (sin(terta 2 f )) 
where X= distance from light to charge, Y=distance from light to grid 
screen, f=geometric factor end Rp=Rp as given in table l. 
This calculation indicated the error was approximately 20% of the 


error shown for large Rp in Figures 18 and 19>. 
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APPENDEX I 
A Computer Data Corporation, Model 1604, Digitsi Computer was 
used to facilitate the numerous end repetitive calculations necessary 
to convert the rew data to useful paremeters. Tha following dis- 
cussion is an explanation of the method end the program utilized to 
accomplish the abevs.o 
Briefly, the programs, as described below, when nsed with the 
library taped BIMED 8 program, will take in the basic data, calculate 
all paremeters used in the thesis, fit the best curves various order 
polynomials to the curve of At versus Ros take the derivative of 
these polynomials aad print out all input and ca.culated parameters. 
(See Table 2) Values of At for given values of a a&é calculated for 
each order polynomial were also printed out to facilitate the plotting 
of these curves. This was necessary to determins whith polynomial 
most closely agreed with the psot of experimental.) determined At and 
Rp 
That portion of the program that was added ts ths BIMED 8 program 
and the reason for doing so is given below. 
Ao Definition of the functions used in the pregrams 
RPRIM(I)e Rp 
TSET(I)= + 
MM = total wumber of deta points used in the program 
PO = atmospheric pressure 
WP s weight of PETN per charge in gramé 
WAZST = waight of lead azide and iead styphnate per charge 
im grams 
TJOULES = energy of explosion per cshargs in joules/gm of PEIN 
VSOUND 2 velocity of sound in standard atmosphere 
WPTOTAL = total equivalents of PETN per charge in grams 
Sena | +3 
(“hu,) 
DELL/ = AA 
ALPHA= 
CONSTI 


TARRIV(I % 
RLAMB(I) = A 
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TRED(I) = 7 

DELTT (1 Sete 

BRATIO(I)=b 

DSCAL(I) = scaled radius 
TSCAL(I) = scaled time 
SLOPE(I)= dat/aR 
SMACH(I)= M P 
APRES(I) =A P in atmospheres 


Be Description of program: 
(1) The following was added directly to the front of the 
BIMED 8 progrem. This portion supplies Ros t% and verious constants 
to the computer and calculates ¢11 perameters not dependent on the 


Slope of the At versus R, curve, 


Pp 
»eJOB*¥LUSK 3 MIN MAX  BIGSHOCK 

PROGRAM BIGSHOK 

DIMENSION RLAMB(50), TARRV(50), DELTIT(50), RPRIM(50), 

1TRED(50), DSCAL(50), TSCAL(50), SMACH(50), APRES(50), 

2SLOPE(50), BRATIO(50), Z(50) 

MM = 36 

READ 810, (RPRIM(I), TSET(I), I = 1,t™) 

READ 701, (Z(I), I= 1,16) 

PO = 14.6959 

WP = .249 

WAZST = .030 

TJOULES = 7690, 

VSOUND = .0343 

WPTOTAL = WP + WAZST*.392/1.73 


SCALE = (WPTOTAL*1.73/2000./453.6 )** 33333 
DELTA = TJOULES*WPTOTAL 
ALPHA = DELTA*100./P0/.68946 )**.33333 


CONSTI = 7./6. 
DO 820 I= 1,MM 
TARRV(I) = RPRIM(1I)/VSOUND 
RLAMB(I) = RPRIM(1I)/ALPHA 
TRED(I) = TSET (I )*VSOUND/ALPHA 
DELTT(I) = TARRV(I) ~- TSET(I) 
BRATIO(I) = RPRIM(1I)/TSET(1)/.001/2.54/12. 
DSCAL(I) = RPRIM(1I)/SCALE/2.54/12. 
820 TSCAL(I) = TSET(1)/SCALE/1000. 
810 FORMAT (£10.4, E20.4) 
701 FORM/T (F5el) 


(2) The input statements of the BIMED 8 were modified 
thusly, to supply Ry and t as its input: 
21 D019 T= 1,NM 
X(I,1) = RPRIM(I)*100.0 
19 X(I,NT) = DELTT(I)*100.0 
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(3) The following statements were edded to the BIMED 8 
immediately after it has determined all the coefficients of the 
polynomiels. This portion of the progrem calculetes the slope of 
the computer determined polynomiels, the Mach number, the overpressure 
end then prints out all parameters for all values of Rp supplied es 
input. It also calculates and prints out the values of at for various 


Rp's as celculated from each polynomiel. This is done, as mentioned 


above, merely for eese of plotting these polynomials. 


PRINT 802 
DO 803 I = 1,16 
Ya 
pO 804 J = 1,N1 
804 Y = Y+B(d )*Z(I )eeJ 
803 PRINT 805, Z(1), Y 
IF (Nl-2) 70,826,826 
826 IF (N1-4) 825,825, 70 
825 DO 850 I=1,.™ 
SLOPE(I) = B(1) + 2.0*B(2)*RPRIM(I) & 3.0*B(3)*RPRIM(I)**2 
\ «L4 4.0#B(4)*RPRIM(I )##3 
SMACH(I) = 1.0/(1.0-VSOUND*SLOPE(I )) 
850 APRES(I) = (SMACH(I)*#2 = 1.0)*CONST1 
PRINT 870, WPTOTAL, DELTA, ALPHA, SCALE 
PRINT 871, Nl 
PRINT 872, B(1), B(2), B(3), B(4) 
PRINT 860 
PRINT 865, (RPRIM(I), TSET(I), TARRV(I), DELTT(I), RLAMB(I), 
ITRED(I), DSCAL(I), TSCAL(I), SLOPE(I), SMACH(I), APRES(I), 
2BRATIO(I), I= 1,MM) 
PRINT 866 
802 FORMAT (1HO SHRPRIM 5X SHDELTT/) 
805 FORMAT (2F10.3) 
870 FORMAT (QHIWPTOTAL® F8.6, ilH DELTA= F9.3, 11H ALPHAs 
1F10.6, 118 SCALES F7.5///) 
871 FORMAT (Q9ODATA FOR I2, 18H DEGREE POLYNOMIAL/) 
872 FORMAT (QHOSLOPES F9.5, 7H + 2.0( F9.5, 13H)RPRIM + 3.0( 
1F9.5, 16H)RPRIM**2 + 4.0( F9.5, 9H)RPRIM**3///) 


860 FORMAT (120HO RPRIM TSET TARRV DELTT 
2RLAMB TRED DSCAL TSCAL SLOPE SMACH 
3APRES BRATIO) 


866 FORMAT (1H1) 
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ALPHA= 26.785693 
SCALE= .00787 


| WPTOTAL= DELTA= 1967.484 - 


°255850 


rte ayes em ew ae 


Semmes = - ee “ ee 


DATA FOR & DEGREE POLYNOMIAL 


SLOPE =._30.204UK7 + 2.0( -.95667)RPRIM + 3.0( eCI889)RPRIMGH2 + 


W.Ol =.CCCIS)RPRIMSS3- 


22.780 


“.* OELTT 


664.140 370.240. 
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RPRIM TSET TARRV RLAMB TREO OSCAL 
.000 -0CO -000 2000 2000 - 20CO .000 
8.220 U8.3CO 239.650 ° 191.350 » 307 2062 34.255 
9-280 47.0C0 270.554 203.554 0 346 -086 38.673 
9.750 68.0CO 284.257 216.257 2364 4087 §=4C.631 
10.370  80.3CO 302.332 222.032. .387 = 103432215 
10.510 80.900 306.414 225.514 392 104 43.798 
10.820 83.4CO 315.452 232.052 e404 . PiCr 45.090 
PT auuin0 84.9CO 334.402 249.502 0428 - 109 47.799 
11.690 88.0CO 340.816 252.816 0436 0113 48.716 
11.780 90-.9CO 343.440 252.540 0440 2116 49.091 
12.400 95.0CO 361.516 266.516. 2463 <198 51.675. 
13.000 110.0C0 379.009 269.0¢9 0485, 014) 54.175 
13.100 110.000 381.924 271.924 0489 + 21h) 54.592 
14.000 128.000 408.163 28C.163 0523 -168 ) Seesnem 
14.100 8 131.8C0 411.079 279.679 0526 0168 58.759 
14.800 184.1CO 431.487 287.387 0552 0184 = 16676 
_14.870 147.100 = 433.528 = 286.428 0555 -188 61.968 
15.450 '151.5CO 450.437 298.937 .-“<577 0154 64.385 
~ 15.730 \ 156.0CO0 458.601 302.601 72587 22C0 éSeo2 
15.9CC 162.000 463.557 301.557 6594 2207 66.260 
36.220 16%.0CO 472.886 308.886 2605 2210 67.59% — 
16.400 165.7CO 478.134 3126434 — 0612 0212 68.344 
16.560 171.000 482.799 311.799 2618 0219 69.010 
16.720 - 173.3CO 487.464 314.164 £62k 6222-6 $2677. 
17.430 190.0CO 508.163 318.163 465) 0243 72.636 
17.850 194.200 520.408 326.208 2666 0249 74.386 
‘19.112 211.000 557.201 346.201 713 2270 79.645 
19.250 215.0CO 561-224 346.224 0719 2275 80.221 
19.310 225.000 562.974  337.97% 72) ~288 80.471 
19.840 -242.8CO 578.426 335.626 274) 0311"  .82.679 
21.300 265.000 620.991 355.991 2795 ~339 88.763 
22.750 295.900 663.265 367.365 2849 °° 379 Su 806 
22.780 293.960 0850376 94931 


mrt a 
—y 


a eg, 
ee ver 
~~, ~t 
oe wr 
Amare ~~ 
wa th. 
aL save 
Wi vi ave 
o-_ te Tra, 
7: © mss 
“hr 7 ii 
do) ee ee 
“> i> ceri 
TT = aah 
oe ier, 
LLL ar. 
oe we 
greta at 
= 4 ™ 
aerrmhlc ; —r. 
we et | ek 
1 -,. 
o» 6 =. 
—- ie sai 
gi Pe 5 ihe 
* “At. 
~a-.* wt ve 
-— @t i “3 6 i 
ye et iF 
: ‘ s+. way 
; f e -*4 ivy 
% : - 
ioe ai :. -—?. 
* si om. 
vs 









_" baa 
"ie 
mE £ 
-— — 
A la 
ei 
atl otbel 
me vil 
e6 @@ 
-” =e 
wrereyy 
™ 995 
<— <a_é 
°°) <—-—_é 
ie 
- ig 
a oe 


wld ale 
wT ffs 














omar iel® 
out OF 
om Lat) 

1) ame 
whbhaawe 
ett i 
“= o8) 
teed! 
“ih 
> ee) 
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RPRIM 


.000 
8.220 
9.280 
9.750 

10.370 
10.510 
10.820 
11.470 
11.690 
11.780 
12.400 
3. O00 
13.100 


. 14.CC00 


14.100 
14.800 
14.870 
15.450 
15.730 
15.960 
16.220 
16.400 
16.560 
16.720 
17.430 


17.850 


W212 
19.250 


19.310 
19.840 


21.300 


2-750 


TSET 

sOCOm 
48.3C0 
7 .0CC 
68.0CO 
80.3C0 
80.9CO 
83.40 
84 .9CO 
88-0CO 
90.9CO 
95.0CO 
110.0C0 
110.000 
128.0CO 
131.400 
144-.1C0 
147.100 
151.5C0 
156.0CO 
162.060 
164 .0CO 
165.760 
reo 
173.3C0 
190.000 
194.2C0 
211.060 
21520C0 
225.000 
242.8C0 
265-0CO 
295.9C0 
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